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ABSTRACT 

Forty-two  in-flight  K    — deuteron  reactions  leading  to  the  final- 
state  A   ,    tt   ,    and  proton  are  investigated.     It  is  shown  that  fewer  of 
these  reactions  involved    2  hyperons  in  an  intermediate  state  than  was 
the  case  for  kaons  interacting  at  rest.      Angular  distributions  in  the 
K-d  and  tt    -A     cm  systems  indicate  P-wave  as  well  as  S-wave 
contribution  to  the  reaction    K     +  n  -*   A     +  it   .      Cross  sections  are 
calculated  for  A    '  s  produced  directly  (30  ±6.  5  mb)  and  those  produced 
via  the  intermediate     2  state  (14.5±4.5  mb). 


-3- 


K   -DEUTERON  INTERACTIONS  IN  FLIGHT 
Paul  G.    White,    Jr. 

Lawrence  Radiation  Laboratory 
University  of  California 
Berkeley,    California 

May  15,    1959 


I.     INTRODUCTION 
In  the  spring  of  1958,    the  Alvarez  group  at  Berkeley  exposed 
a  15-inch  bubble  chamber  filled  with  liquid  deuterium  to  a  beam  of 
slow  (about  250-Mev/c)  negative  K  mesons  in  order  to  study  the 
elementary  K-nucleon  reactions.      Of  the  various  types  of  interactions 
observed  in  the  chamber,    one  lent  itself  most  cleanly  to  kinematic 
analysis,    namely  that  in  which  a  kaon  and  a  deuteron  interacted  to  give 
a  negative  pion,    a  proton,    and  a   A     hyperon  which  subsequently  de- 
cayed in  the  chamber  into  a  negative  pion  and  a  proton.      This  is  the 
only  K-deuteron  reaction  (except  for  elastic  K-d  scatterings)  in  which 
all  the  final-state  particles  leave  visible  tracks,    which  makes  the 
solution  for  the  kinematic  variables  overdetermined  by  the  four  con- 
straints of  energy  and  momentum  balance.      This  type  of  event  was 
given  the  numerical  designation  769.      Upon  analysis,    a  small  fraction 
of  these  events    (about  15%)  were  found  to  have  occurred  while  the 
kaon  was  in  flight.     It  is  this  subclass  of  events  which  is  the  subject 
of  this  paper. 
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II.      KINEMATIC  ANALYSIS 

A.  Hand  Analysis 

Included  in  the  general  category  of  Type  769  events  are  four 
subclasses,    all  of  which  have  the  same  general  appearance  in  the 
chamber  (see  Fig.    1):     A     produced  at  rest,    A     produced  in  flight, 

2     produced  at  rest,    and  S     produced  in  flight.      A  S     production 

0  0  -19 

appears  similar  to  a  A     because  the    S     decays  in  ~  10  second  into 

a     A     and  a  gamma  of  about  7  5  Mev  energy.     At  this  decay  rate,    even 

0  -8 

if  the    S     had  the  velocity  of  light,    it  could  move  only  ~  10        centimeter 

before  it  decayed  into  a   A     and  a  gamma.      The  only  outward  difference 
between  the  two  types  of  events,    then,    is  that  the    S     production  has 
an  energy  and  momentum  unbalance  of  about  7  5  Mev,    due  to  the  missing 
gamma.     All  Type  769  events  were  analyzed  by  a  hand-analysis  tech- 
nique which  took  as  inputs  the  azimuthal  and  dip  angles  of  each  track, 
the  momenta  as  determined  by  track  curvature,    or- -for  stopping  tracks 
the  momentum  as  determined  by  track  length.      These  data  were  ob- 
tained as  output  from  an  IBM  650  program  which  took  digitized  data 
from  two  stereo  photographs  of  the  bubble  chamber  event  and  printed 
out  the  azimuthal  angle,    dip  angle,    track  length,    momentum  (from 
track  curvature),    and  included  angle  between  each  pair  of  tracks  in 
the  event.      The  parameters  were  then  adjusted  by  trial  and  error 
until  a  consistent  solution  was  found  which  balanced  momentum  and 
energy  at  the    A    -production  and  decay  vertices.      Usually,    this 
technique  would  distinguish  unambiguously  between  the  four  possible 
subclasses. 

B.  Computer  Analysis 

Although  it  was  found  possible  to  separate  out  the  A    '  s  pro- 
duced in  flight  by  this  method,    there  are  so  many  degrees  of  freedom 
in  this  type  of  event  that  it  is  difficult,    if  not  impossible,    to  find  by 
hand  analysis  a  "best"  solution  which  balances  both  momentum  and 
energy.      Fortunately,    an  IBM  704-computer  program  designed  to 
perform  this  type  of  analysis  has  been  developed  by  Dr.    Horace  Taft 
and  Peter  Berge.        The  events  involving  in-flight  A     production  were 
analyzed  with  this  program  in  order  to  determine  the  kinematic 
parameters  more  accurately  than  was  possible  by  hand  analysis. 
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Fig.    1.      Representative  Type  769  event  in  flight. 
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Input  variables  to  the  program  for  each  of  the  four  tracks  at 
the  production  vertex  were  (a)  the  azimuthal  angle;  (b)  tangent  of  the 
dip  angle;  and  (c)  l/p  cos  \,  where  X  is  the  dip  angle  and  p  is  the 
track  momentum  as  determined  by  curvature,  by  range  (for  stopping 
tracks),  or--in  the  case  of  the  neutral  A  --by  fitting  the  parameters 
at  the  A  decay  vertex  on  a  kinematics  chart.  These  variables  were 
chosen  because  their  errors  have  a  nearly  Gaussian  distribution. 
The  variance  of  each  of  the  twelve  variables  was  also  an  input  to  the 
program. 

The  problem  confronting  the  computer  is  to  calculate  a  fit  to 

2 
the  input  parameters  which  minimizes  the    x    --subject  to  the  side 

conditions  that  the  x,    y,    and  z  components  of  momentum  must  be 
balanced,    and  that  energy  must  be  conserved.      A  geometrical  in- 
terpretation of  the  problem  would  be  that  in  the   12-dimensional  space 
representing  all  the  possible  values  of  the   12  input  parameters,    there 
is  an  8-dimensional  "surface"  which  includes  only  the  values  of  the 
parameters  allowed  under  conservation  of  momentum  and  energy. 
The  input  variables  locate  a  point  in  the   12-dimensional  space  which 
will  not  in  general  be  located  on  the  surface  of  allowed  situations, 

and  this  point,    together  with  the  input  variance  for  each  variable, 

2 
determines  the  value  of  x      f°r  each  point  in  the   12-dimensional  space, 

The  problem  then  is  to  locate  the  point  on  the  8-dimensional  surface 

at  which  x      is  a  minimum.      The  problem  is  solved  by  the  introduction 

of  Lagrange  multipliers  for  the  four  constraint  equations,    giving   16 

equations  in  16  unknowns.      Because  of  the  nonlmearity  of  the  constraint 

equations,    an  iterative  procedure  is  used  which  assumes  moderate 

linearity  in  the  region  containing  the  points  corresponding  to  the 

initial  and  final  values  of  the   12  parameters.      The  machine  will 

2  .  .  2  . 

continue  to  take  "steps"  to  minimize  x      until  the  gradient  of  x      is 

less  than  a  given  maximum  acceptable  value  and  the  total  momentum 

and  energy  unbalance  is  less  than  1  Mev.     At  this  point,    the  computer 

prints  out  the   12  fitted  variables,    a  new  variance  for  each,  and  the 

value  of  x      f°r  the  event.      Since  in  most  events  four  constraints  were 

applied  to  the  problem,    an  average  x      of  4  was  expected.      In  some 
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cases,    however,    the  proton  from  the  production  vertex  had  such  low 
momentum  that  its  track  length  was  too  short  to  be  visible.      The 
computer  was  then  forced  to  use  up  three  of  the  constraint  equations 
in  calculating  the  missing  momentum  vector,    leaving  the  problem 
overdetermined  by  only  one  degree.      For  these  cases,    an  average 
X      of  1  was  expected.      The  x      f°r  four -constraint  cases  averaged 
4.44;  for  the  one -constraint  cases  the  average  was   1.43. 

Out  of  47  events  that  were  analyzed  by  the  program  as  being 

probable  A     events  in  flight,    the  analysis  indicated  that  three  actually 

Z 
occurred  at  rest.      Two  others  with  anomalously  large  x       s,    it  was 

concluded,    were  actually   2     events. 

Standard  deviations  in  calculated  momenta  were  typically 

±  5  Mev/c  for  four -constraint  cases,    and  ±20  Mev/c  for  one -constraint 

cases.     Angles  were  determined  to  within  about  l/2  degree  in  all 

cases  except  for  the  calculated  azimuth  and  dip  of  invisible  protons, 

which  typically  had  uncertainties  of  about  20    . 


III.      RESULTS 

A.      Direct  vs.    Indirect  A     Production 
1.     At-Rest  Case 

Analysis  of  the  at-rest  Type  769  events  indicates  that  a  large 
fraction    (about  60%)  of  the    A     hyperons  were  produced  by  a  two-step 
process  in  which  the  kaon  reacts  with  one  of  the  nucleons  to  give  a 
2     or  2     and  a    it.        The     2    hyperon  then  interacts  with  the  other 
nucleon  to  give  a   A     and  a  proton.      The  evidence  for  this  process  is 
that  the  momentum  distribution  of  the  pion  from  the  A     production 
vertex  shows  two  pronounced  peaks,    one  at  about  250  Mev/c,    and  a 
larger  peak  at  about  185  Mev/c   (Fig.    2).      The  peak  at  250  Mev/c  is 
interpreted  as  being  due  to  those  events  in  which  the  A     was  pro- 
duced directly  from  a  K  -neutron  interaction,    the  proton  merely 

being  a  "spectator"  in  the  process.      It  is  similar  in  form  to  the 

3 
distribution  calculated  theoretically  by  Fujii  and  Marshak.        The 

second  peak  is  considered  to  have  arisen  from  the  two-step  process 

0  +  4 

involving  a     2     or  2     in  the  intermediate  state,        since  it  corresponds 
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Fig.    2.      Momentum  distribution  for  pions  for  228  at-rest 
events. 
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to  the  calculated  pion -momentum  distribution  for  a    2    and  a    tt 
produced  by  a  kaon  in  deuterium,    yet  the  over-all  kinematics  for  each 
event  fits   A     production,.      The  at-rest  Type  769  events  were  thus 
divided  into  two  categories  on  the  basis  of  the  momentum  of  the  pro- 
duction pion.      Those  with  pion  momenta  greater  than  215  Mev/c  are 
considered  to  involve  direct  A     production  and  those  with  pion 
momenta  less  than  this  are  considered  to  involve  the  two-step  or  in- 
direct production.      The  indirect  at-rest  events  taken  on  the  basis  of 
this  separation  outnumber  the  direct  by  a  ratio  of  about  8/5. 

2.  In-Flight  Problem 

In  order  to  determine  this  ratio  for  the  in-flight  events,    clearly 
a  new  criterion  must  be  found  for  separating  direct  from  indirect 
events.      The  in-flight  events  differ  from  the  at-rest  case  in  that  the 
center   of  mass  is  moving  in  the  laboratory  system  of  reference  and 
the  kinetic  energy  of  the  kaon  brings  a  variable  amount  of  additional 
energy  into  the  center-of-mass  system,    some  of  which  will  go  into 
increasing  the  pion  momentum  above  what  it  would  have  been  had  the 
event  occurred  at  rest. 

3.  Method  of  Separation 

If  one  assumes  that  the  impulse  approximation  is  correct  (l.  e.  , 
in  the  direct  events  the  proton  is  merely  a  "spectator"  in  the  interaction^ 
the  distribution  of  proton  momentum  in  the  laboratory  system  can  be 
calculated  (see  Appendix).      This  predicted  momentum  distribution  is 
compared  with  the  experimental  distribution  in  Fig.    3.      The  comparison 
shows  that  there  are  too  many  protons  with  high  laboratory-system 
momenta.      One  obvious  interpretation  of  this  fact  is  that  these  high- 
momentum  protons  were  the  ones  involved  in  the     S-to- A  conversion 
process.      Since  this  process  is  an  exothermic  reaction  between    2!  s 
and  nucleons,    the  protons  resulting  are  expected  to  have  fairly  high 
momenta.      One  might  establish  the  criterion  that  all  events  with 
production  proton  momenta  less  than  200  Mev/c  were  direct  events 
while  all  others  were  indirect.     Since  about    20%  of  the  direct  events 
that  occurred  at  rest  had  proton  momenta  greater  than  200  Mev/c  (see 
Fig.    4),   the  classification  of  events  on  this  basis  establishes  the  lower 
limit  to  the  number  of  direct  events. 
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Fig.    3.      Comparison  of  predicted  and  experimental 
momentum  distributions. 
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Fig.    4.      Momentum  distribution  for  protons  for  97  direct 
events  (at  rest). 
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4.      Results 

The  separation  of  events  by  proton  momenta  indicates  25 
direct  events  and  17  indirect  events.      It  is  clear  that  the  ratio  of  in- 
direct  A    '  s  to  direct  A    '  s  when  the  K  interacts  in  flight  is  consider- 
ably less  than  the  ratio  when  the  K  interacts  at  rest. 

B.      Angular  Distributions 
1.      K-Deuteron  Center-of-Mass  System 

If  one  regards  the  process  K  +  d -»•  tt     +  A     +pas  one  in- 
volving a  two -body  initial  state  and  a  three -body  final  state,    one  can 
ask  what  partial  waves  in  the  K-d  system  are  involved  in  the  reaction. 
Since  the  momenta  of  the  incoming  kaons  is  low    (  <  250  Mev/c),    there 
is  little  possibility  that  angular  momenta  higher  than  P-wave  are  in- 
volved.      Under  this  hypothesis  the  angular  distribution  of  any  of  the 
reaction  products  must  have  the  form 

42-   =   1  +  AcosQ  4  Bcos20, 
aki 

where     G    is  the  angle  between  the  incoming  kaon  and  the  outgoing 

particle  in  the  K-d  center-of-mass  system.      If  only  S-wave  angular 

momentum  is  involved,    then  A  =  B  =  0  and  the  angular  distribution  of 

all  reaction  products  must  be  isotropic  with  respect  to  the  kaon.      The 

histogram  of  cos  9    ,,(K-d,    c.  m.  )  for  25  direct  events  chosen  on  the 

TXlS. 

basis  of  a  laboratory-system  proton  momentum  less  than  200  Mev/c 

is  shown  in  Fig.    5.      A  least-squares  fit  of  these  data  to 

1  +  A  cos    6    ^  B  cos    0    „      gives    A  =   1.04±.52,      B  =  0.68±.85. 

TTiS.  7TXV 

Values  of  A  and  B  for  all    events,    direct  and  indirect,    in  the  K-d 
system  are 

A  =  0.58±.31, 

B  =  0.33±.51. 
Nonzero  values  of  A  and(or)  B  for  the  distribution  in  the  K-d  center- 
of-mass  system  does  not  necessarily  imply  the  existence  of  P-wave 
angular  momentum  in  the  K-neutron  interaction,    however.      It  has 

been  calculated  that  200-Mev/c  kaons  interacting  only  by  S  channel 

5 
with  neutrons  would  give  a  value  A  =  0,17.        This  may  be  seen 

aualitativelv  bv  considering  that  since  the  proton  must  have  a  low 
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Fig.    5.  Distribution  of  cos.0    K  in  K-d  center-of-mass 
system  for  25  direct  events    (P     <  200  Mev/c). 

ir 


14- 


laboratory-system  momentum,    it  tends  to  be  directed  backward  in 
the  K-d  center-of-mass  system;  this  backward  proton  momentum 
must  be  balanced  by  having  the  pion  and  A     tend  slightly  forward  with 
respect  to  the  kaon. 
2.    tt-A     Center-of-Mass  System 

In  order  to  determine  whether  there  is  a  P-wave  component  in 
the  K-neutron  interactions    one  would  like  to  examine  the  distribution 
of  the  tt-K  angle  in  the  K-neutron  center-of-mass  system  instead  of 
the  K-deuteron  system,      Since  the  neutron  is  bound,    it  is  in  an  un- 
physical  region  and  the  K-neutron  system  cannot  be  measured  directly. 
To  the  extent  that  the  neutron  rest-mass  energy  is  large  compared  with 
the  deuteron  binding  energy,    the  tt-A     center-of-mass  system  should 
have  nearly  the  same  velocity  and  total  center-of-mass  energy  as  the 
K-neutron  system.      In  the  impulse  approximation,    then,    the  distribution 
of  the  tt-K  angle  in  the  tt-A     system  should  be  the  same  as  that  distri- 
bution in  the  K-neutron  system.      The  effect  of  the  inaccuracy  of  this 
approximation  should  be  to  smear  the  angular  distribution  out  some- 
what from  what  it  would  have  been  in  the  K-neutron  system,    but  it 
should  not  introduce  any    "artificial"  asymmetries.      The  distribution 
for  the  25  direct  events  is  shown  in  Fig.    6.      The  parameters  calculated 
from  this  distribution  are 

A  =  0.7  5  ±.44, 
B  =  0.43±.77. 
From  this,    it  appears  that  the  direct  A     production  proceeds  partially 
by  P-wave  even  at  these  low  kaon  momenta. 

Co      Cross  Sections 
In  estimating  cross  sections  for  the  direct  and  indirect  pro- 
duction in  flight,    it  is  necessary  to  exclude  events  that  occur  near  the 
chamber  boundaries,    since  the  detection  efficiency  for  those  events 
is  probably  low  and  its  value  is  unknown.      Assuming   100%  detection 
efficiency  for  those  events  occurring  near  the  center  (in  the  "defined" 
chamber  region),    it  is  then  necessary  to  calculate  the  total  kaon  path 
for  all  kaons  having  tracks  that  enter  the  defined  chamber.      Since  the 
tracks  of  kaons  that  passed  completely  through  the  chamber  have  not. 
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been  analyzed,    it  is  necessary  to  exclude  a  second  class  of  events 
from  consideration  in  determining  cross  sections:     those  which 
occurred  in  such  a  chamber  location  and  with  such  kaon  momenta 
that,    had  the  kaon  not  interacted  in  flight,    the  kaon  would  not  have 
stopped  within  the  chamber. 

The  total  kaon  path  length  in  various  momentum  intervals  for 
the  class  of  kaons  which  either  stopped  within  the  chamber  or  would 
have  stopped  had  they  not  interacted  in  flight  has  been  calculated  in 

connection  with  another  experiment.        The  average  cross  section 

2 
(in  cm    )  for  a  given  reaction,    in  a  given  momentum  interval   p  +  Ap 


and  p,    is  obtained  by 


O      -  N  /pLp, 
P  P    V 


where  N     is  the  number  of  events  that  occurred  when  the  kaon 

P 
momentum  was  in  the  range     Ap,    p  is  the  number  of  deuterons  per 

cubic  centimeter,    and  L     is  the  total  path  length  in  centimeters 

traveled  by  kaons  with  momenta  between  p+  Ap  and  p.      Cross  sections 

for  direct  and  indirect  Type  769  events  are  listed  in  Table  I. 


Table  I 

Cross  sections  for  direct  and  indirect  Type  769  events' 

Direct  Indirect 


K  momentum 

No.    of 

Cross  sec 

tion 

No. 

of 

Cross  section 

(Mev/c) 

events 

(mb) 

eve 
3 

nts 

(mb) 

50  -    100 

3 

49  ±28 

49  ±28 

100  -   150 

5 

21±    9 

3 

13±    8 

150  -  200 

13 

Av. 

25±    8 

5 

4 
10 

Av 

10±    5 

50  -    200 

21 

30±    6. 

14.5±4.5 

3. 

Errors  are  statistical  only. 
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IV.      CONCLUSIONS 

In  interpreting  the  results  of  this  experiment,    it  is  necessary 
to  know  whether  the  change  in  the  ratio  of  direct  to  indirect   A 
production  is  due  to  a  decrease  in  the  percentage  of    2  hyperons  that 
convert  to   A    '  s  or  to  an  increase  in  the  ratio  of  A    '  s  produced 
directly  relative  to  the  number  of   2  hyperons  produced. 

On  the  basis  of  the  models  of  the  conversion  reaction  considered 

4 
by  Karplus  and  Rodberg,      one  would  not  expect  the  conversion  process 

7 
to  be  materially  changed  by  allowing  the  kaons  to  interact  in  flight. 

There  is  evidence  that  in  K    -p  reactions  in  hydrogen,    direct   A 

production  is  substantially  increased  relative  to    2  production  when 

o 

the  kaon  interacts  in  flight  (300  to  400  Mev/c)  instead  of  at  rest. 
The  most  likely  explanation  for  the  change  in  the  direct -indirect  ratio 
then  appears  to  be  that  the  K    -n  reaction  also  is  more  likely  to 
produce     A    '  s  directly  when  the  kaon  interacts  in  flight. 
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APPENDIX 

In  the  calculation  of  the  matrix  element  for  direct  A 
production,   the  following  assumptions  are  made: 

(a)  The  kaon  and  neutron  interact  at  a  point. 

(b)  The  pion  and  A     are  created  at  the  same  point. 

(c)  No  secondary  interactions  involving  the  proton  occur 

(i.  e.  ,    no  TT-p  or  A-p  scattering). 

The  variables  are  expressed  in  units  such  that     fi  =  c  =   1.      The  origin 

of  all  spatial  coordinates  is  the  deuteron  center  of  mass,    and  all 

momentum  vectors  are  with  respect  to  the  laboratory  system  of 

reference,    so  that  the  conditions  r     =r     and  P,   =  P    +P    +PA    apply. 

pn  k        p         tt        A 

The  initial  state  may  then  be  described  by 


ifj  (r    ,  r    ,  r,  ) 
xi     p      n      k' 


<2ir) 


6(P-Pk)    exp(iP°  ?k)4>d  (P1)  eV  '   (rp-rn)d3Pd3P', 

(1) 


where    cj>  ,     is  the  momentum  representation  of  the  deuteron  wave 
function.      The  final  state  is  given  by 


*f  ^A'V  rp»  =    - 


^72 jff  6<Q1-PA»  6(Q2-P,»  6<Q3-Pp» 


exp 


i  (Qj  •    rA  +  Q2  •    r^  +  Q3  •    r    ) 


3         3         3 
dJQ1d°Q2d°Q3 


(2) 


The  matrix  element  for  direct  A     production  is  given  by 

/  ij;/  T  i|i.    dT,  (3) 

V  x 

where     T    is  the  interaction  operator  specified  in  this  approximation  by 

T  =  6(rk-rn)6(vrA)6(rA-rn).  (4) 


-20- 
Then,    substituting  Eqs.    (1)  and  (2)  into  (3),    one  has 
M  =  (-^T5/2  /v  S(Q1-PA)8(Q2-Pu)6(Q3-Pp)exp[Q1-   rA+  Q2  •   r^C^  rp)]  T 


6(P-Pk)c|)d(P,)exp[i(P-  rk+P-  .  r     -P'  •  r    )]  d3Pd3P'd3Qj 


3  3  3         3         3         3         3 

d°Q9dJQ_d    r,  d    r    d    r    d    r    d    rA.  (5) 

2  3         k        n        p         it        A 

Integrating  over  P,    Q,,    Q    ,    Q~  gives,    from  the  delta  functions, 


M  ■  (7^3/2     }v  «*P[-»(PA-  *A+P»-  -W  V]T  V1"' 


exp  [i  (P    *   r,   +  P1  (r     -r    )1  ■    d3P1d3r1  d3r    d3r    d3r    d3rA. 
^Lkk  pn/J  knp-rrA 

(6) 
Substituting  for  T  and  integrating  over  r,  ,    r    ,    and  r.    gives 


M  =  yr,      /        expfi  (P    •   r    +  P'  •   r     -  P'  •   r     -  P.  •   r     -  P    '   r    -  P   •   r    )1 

(2tt)3/2     Jv  ^LXkn  p  n         A       n         tt       n        p       p'1 


A,   (P'  )d3P'd3r    d3r    ,  (7) 

Yd  n         p 


M=7-V"2    L    exp[i{(pk-pA-Vp'»'r„+<p'-pp»'rP}1*d(p,) 

(2tt)  Jv 


d3P'd3r    d3r      .  (8) 

n         p 


2.1 


From  the  side  conditions  r      =  r     and  P,=P.+P     +P,    one  obtains 

n         p  k  A  n  p 


M  = ^-r,      /  exp    [  2i    (P'  -  P    )  •  ■  r     1  <j>,   (P1  )  d3r   d3P'  . 

(2u)3/2     ^V  P  P         d  P 


(9) 


Integrating  over  r     gives 


M  =  —^-r?      \         MP'  -  P   )  4>H  (P'  )d3P'    =  4>H(P   >.  (10) 

(2tt)3/2     Jv  P       d  d      P 

The  momentum  distribution  for  protons  in  the  laboratory  system  for 
direct  events  is  given  by 


dN  -,     dp  _     dp 


=  |M|"   —      =   l<MPJI       TO  .  (11) 


dP  ■       '       dP  IYdv    p"       dP 

P  P  P 


where   -rp         is  the  phase  space  factor  describing  the  density  of  states, 

P 
This  distribution  is  plotted  in  Fig.    3„ 
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